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Electron-hole (e-h) capture in luminescent conjugated polymers (LCPs) is modeled by the dissipa- 
tive dynamics of a multilevel electronic system coupled to a phonon bath. Electroinjectcd c-h pairs 
are simulated by a mixed quantum state, which relaxes via phonon-driven internal conversions 
to low-lying charge-transfer (CT) and excitonic (XT) states. The underlying two-band polymer 
model reflects PPV and spans monoexcited configuration interaction singlets (S) and triplets (T), 
coupled to Franck-Condon active C=C stretches and ring-torsions. Focusing entirely upon long 
PPV chains, we consider the recombination kinetics of an initially separated CT pair. Our model 
calculations indicated that S and T recombination proceeds according to a branched, two-step 
mechanism dictated by near e-h symmetry. The initial relaxation occurs rapidly with nearly half 
of the population going into excitons (Sxt or Txt)> while the remaining portion remains locked 
in metastable CT states. While formation rates of Sqt an <i T^CT are nearly equal, Sxt 1s formed 
about twice as fast Txt m concurrence with experimental observations of these systems. Further- 
more, breaking e-h symmetry suppresses the XT to CT branching ratio for triplets and opens a 
slow CT— » XT conversion channel exclusively for singlets due to dipole-dipole interactions between 
geminate and non-geminate configurations. Finally, our calculations yield a remarkable linear re- 
lation between chain length and singlet /triplet branching ratio which can be explained in terms of 
the binding energies of the respective final excitonic states and the scaling of singlet-triplet energy 
gap with chain length. 



I. INTRODUCTION 

Since the discovery of electroluminescence (EL) in poly(p-phenylene vinylene) PPV,(0) considerable efforts have 
gone in development and commercialization of light-emitting diodes (LEDs), based on luminescent conjugated poly- 
mers (LCPs). While much progress has been made in improving emission efficiency and color control, the detailed 
physics of charge transport and light generation in these materials is relatively poorly understood in comparison with 
inorganic semiconductors. The one-dimensional derealization of electrons and substantial electron-hole and electron- 
lattice interactions render the delineation of neutral and charged excitations in LCPs very complex and demanding. 
Moreover, a dynamical description is needed in order to capture the transient relaxation processes that occur when 
charged species approach each other, recombine and decay. 

Photogenerated singlet e-h pairs in LCPs typically relax to highly emissive S\ excitons prior to decay. In contrast, 
electroinjected electrons (e - ) and holes (h + ) are not spin correlated and can combine to form both singlet and triplet 
excitons. 
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Because the decay of the triplet exciton T\ is nearly radiationless and relaxation to S occurs via spin-orbit coupling, 
quantum efficiency in electroluminescence (tjel) is only a fraction of that for photoluminescence (tjpl). If the rate of 
e-h capture is spin- independent, S± and T\ excitons are formed in a 1:3 ratio as dictated by spin degeneracy. In such a 
case, 75% of the bound e-h pairs are non-emissive leading to a theoretical maximum of tjel = 0.25 r\ph- Nonetheless, 
efficiencies of up to rjpr = 0.50 rjp l have been achieved independently on PPV-based LEDs in the laboratories of 
Heeger|0) and Friend(|3). From this it has been inferred that in organic light-emitting polymers, singlet e — h capture 
is intrinsically more efficient than the respective triplet process (ks > kp). 

If we assume that the formation kinetics is first order such that the singlet and triplet exciton population formation 
rates are proportional to the total population of charge-transfer (CT) states in the system, then the formation cross- 
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section a s fov singlet excitons as given by the ratio of singlets to total excitons by 

N s (t) r 



(2) 



N s (t) + 3N T (t) r + 3 

where r — ks/kr = os /&t- In order to obtain the nearly 64% EL efficiencies reported by various groups|2 0) in 
various organic polymer LED devices, r must be in the range of 3 < r < 5. Such ratios were also obtained for a wide 
range of organic LCPs via photoabsorption/detection of magnetic resonance (PADMR) experiments by Wohlgenannt, 
at al 10 0) which measures the population of photo-generated intrachain parallel and anti-parrallel polaron pairs in 
the presence of a magnetic field. This later work|0) is particularly significant in that it establishes a universal linear 
scaling relationship between r and inverse conjugation length, which holds for a wide range of LCP systems. 

Various mechanisms favoring the formation of Si versus T\ have been proposed for both interchain and intrachain 
c-h collisions. Using Fermi's golden rule, Shuai, Bredas et al.|2@) indicate that the S cross section for interchain 
recombination can be higher than the triplet one due to bond-charge correlations. Wohlgenannt et aL0) resort to a 
similar model of two parallel polyene chains. Both of these works ne glec t vibronic and relaxation effects. In simulating 
the intrachain collision of opposite polarons, Kobrak and Bittner ifflulft) show that formation of S% is enhanced by the 
near-resonance with the free e-h pair. The result reflects the fact that spin-exchange renders T\ more tightly bound 
than Si .ijllt ll^) and hence more electronic energy must be dissipated by the phonons in the formation of the former. 
The energy-conservation constraints in spin-dependent e-h recombination have been analyzed by Burin and Ratner l|l3l) 
in an essential-state model. The authors point out that nonradiative processes (internal conversion, intersystem 
crossing) must entail C=C stretching vibrons since these modes couple most strongly to 7r — > 7r* excitations. Recent 
work by Tandon et al. suggest that irrespective of the recombination process, interchain or intrachain, the direct 
transition to form singlets should always be easier than triplets due to its smaller binding energy relative to the 
triplet. l)l4h . Lastly, a comprehensive review of this work detailing the experiments and summarizing the theory of 
this effect is forthcoming. 1|15() 

In this work, we describe e-h capture in LCPs in terms of the dissipative dynamics of a multi-level electronic system 
coupled to a phonon bath. Our model is based upon a two-band polymer model introduced by Soos, Mazumdar, et 
al l)l6l : Il7j) augmented by the coupling of electronic excitations to a bath of vibrational modes. First, we review a 
time-dependent model for simulating via quantum chemical methods the transient relaxation dynamics of an excited 
state in an extended system. l(l8t I20I) The methodology we develop herein is applicable to both excitonic transfer (i.e. 
Forster) as well as charge-transfer states in a general conjugated polymer system. At the present, we restrict our 
attention to mono-excited closed-shell systems; however, the approach can be extended to the general case for radicals. 
Secondly, we examine a specific physical process involving the collision and relaxation of injected electron-hole pairs in 
luminescent conjugated polymers, such as PPV. ljl9f) Our calculations underscore the role that vibronic coupling and 
electron-hole symmetry play in determining the singlet-triplet branching ratio as the CT state relaxes to form bound 
electron/hole excitonic pairs. Finally, we explore the relation between chain length and the singlct/triplc branching 
ratio. We show that our methodology reproduces the remarkably universal linear relation between conjugation length 
and r as evident in a wide range of conjugated polymer materials. To our knowledge, ours is the first molecular based 
model which accurately predicts this linear relation for a specific molecular system and provides a rationalization for 
this trend based upon the variation of the exchange energy with increasing chain length. 



II. METHODOLOGY 

A. Two-band polymer model + phonons 

The Hamiltonian for the coupled system is 

H = H e i + H p h + H e i^ p h 

mn 

fjLU 

The electronic Hamiltonian H e [ represents the configuration interaction of localized singlet of triplet configurations 
|m) = |mm) with a valance hole in repeat unit rn and a conduction electron in repeat unit m. For example, the \mm) 
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geminate monoexcited configuration for a 6 repeat-unit system is 

\mm) = (0-0-0-0-0-0) (4) 

where as 

\mn) = (0 - - - - - 0) (5) 

represents a charge-transfer configuration. Our notation is such that (— — ) represents a neutral site in the ground- 
state configuration, (— — ) denotes a site with a hole, (— — ) a negatively charged site with an excess electron, and 
(— — ) a neutral excitonic site. 

System specific information is incorporated by using valance and conduction band Wannier functions (WF) |m) 
and \m) as the single particle basis. For the single-particle band-structure we use a Hartree-Forck Pariser-Parr-Poplc 
(PPP) approximation which can be parameterized easily for a variety of luminescent conjugated polymer systems, and 
construct maximally optimized Wannier functions. The single-particle terms, F^ n , are derived at the ground state 
equilibrium configuration, , = 0, from the Fourier components f° and f r of the band energies in pseudomomentum 
space. 

F mn = S—i{ m \f°\ n ) - 5mn(rn\f°\n) 

— ^mnfm-n ~ ^mnfm-n (6) 

For extended 7r-systems with electron- hole symmetry, the localized electron/hole energy levels and transfer parameters 
are related by f r = —f r and reflect the cosine-shaped valance and conduction bands of half- width f° and centered at 
±/° respectively. 

The two-particle interactions are spin dependent with 

Vmn = -{mri\\nm) (7) 
V£ n = VT n + 2{mn\\mn). (8) 

with 

(mn\\ij) = JdlJ d20 m (l)*<^(2Ml2)<Ml%(2) 

With the exception of geminate WFs, orbital overlap is small such that the two-body interactions are limited to 
Coulomb and exchange integrals reflecting e-h attraction and spin-exchange coupling non-geminate configurations 
such as 

J\m-n\ = (mn\\mn) = J a (l - r/ro)^ 1 (9) 
K\m—n\ = (mn\\nm) = K exp(—r/r ) (10) 

and transition dipole-dipole integrals coupling only geminate singlet electron- hole pairs, 

D\ m - n] = D (r/r ) 3 . (11) 

In the absence of e-h symmetry, dipole-dipole coupling between singlet geminate and non-geminate configurations 
is introduced through this term. As we shall explore later in this paper, this opens a channel for coupling singlet 
charge-transfer states to singlet excitonic states but does not allow for coupling between the triplet CT to XT states. 

Finally, the electron-phonon coupling term, H e i- p h describes the local variation of the single-particle band gap and is 
taken to be linear a set of phonon normal-mode coordinates which are localized within a single unit-cell. The linear 
coupling strength for both the electron and holes, (df jdq^) = —(df/dq^o, is adjusted empirically as to reproduce 
the the vibronic progressions observed in the single-photon absorption and emission spectra. Correspondingly, the 
vibrational term, H p h models weakly dispersec optical phonon branches in the frequency ranges of the dominant 
Franck-Condon active modes observed in the experimental spectra. For the case of PPV considered here, the phonon 
term consists of two sets of local harmonic oscillators with weak nearest-neighbor coupling, representing two dispersed 
optical phonon branches, centered at 1600 and 100cm -1 , respectively. These frequencies roughly correspond to C=C 
bond stretches and ring-torsional motions (librations), which dominate the Franck-Condon activity in PPV and related 
polymers. |U |3 E3) Huang-Rhys factors (dimensionless el-ph coupling constants) for Si of the isolated monomer 
were subsequently set to 0.8 for the high-frequency mode and 10 for the low- frequency mode. 
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B. Diabatic Representation 

Having set up the primitive Hamiltonian and parameterized it according to a model physical system, we can obtain 
the diabatic electronic/ vibrational structure by separately diagonalizing H e i and H p h yielding a series of vertical 
excited states with energies e° (from = 0) and a normal modes with frequencies u x i. Written in the diabatic 
representation 

H = 5»><a| 



+ E^Q«|a>(b|. (12) 

The diabatic coupling term g° b ^ is simply the original clcctron-phonon coupling tensor transformed into this new 
representation. The density of states (DOS) for vertical singlet and triplet excitations from the ground state for a 
PPV model with 32 repeat units (PPV32) is shown in Fig.^ a s parameterized from the 7r-band structure and Wannier 
functions of the extended polymer. Although the Stokes shifts from adiabatic relaxation of the excited states is missing 
in this calculation, the lowest singlet Si and triplet T\ excitons differ significantly in terms of their binding energies. 
Here, Ti is roughly lcV lower than the Si in general agreement with experimental and theoretical estimates. We 
can also use the diabatic representation to calculate the relaxed excited states and their respective potential energy 
curves. This we achieve by variational minimization of the adiabatic energy e a = (a|/f |a) of each individual state 
according to 

de 

floof + w|Q c = 0. (13) 

The relaxed excited states are eigenstatcs of H e i + H p h at state specific conformations, Q^, a 0; displaced from the 
ground state at = 0. Such displacements modify the o ptic al spectra by moving states into and out of regions of 



strong Franck-Condon coupling to the ground state. I|26t l27j) Thus, we can test the parameterization of the model 
by computing the Condon spectral density for single photon absorption. The cumulative spectral density of PPV32 
is given in Fig. El for both the vertical and Condon approximation. Individual electron-vibrational transitions are 
taken to give Lorentzian lincshapes of half-width O.OleV. The oscillator strengths for the S — ► S a vertical transitions 
were computed by assuming that only geminate and nearest neighbor e-h pairs are dipole coupled to S . Since Si is 
approximately and all-symmetric combination of such configurations, the bulk of the spectra density is concentrated in 
this state. Furthermore, the S a — > Si vibronic band is peaked in the — > 1 C=C stretching feature in agreement with 
the experimental spectra. Vibronic structure from librations is smeared and produces spectral broadening proportional 
to Huang-Rhys factors times phonon frequencies. The low-frequency coupling also determines the « 0.15eV Stokes 
shift of the 0^0 C=C stretch feature from the adiabatic origin of Si at 2.34cV. 



C. Relaxation dynamics 

Separated CT states are not eigenstates of the diabatic Hamiltonian and will evolve in time according to the time- 
dependent Schrodinger equation. Unfortunately, a complete description of the vibronic dynamics is unfeasible due 
to the enormous size of the state-space. However, we can consider the phonons as a dissipative finite-temperature 
bath for the electrons and derive reduced equations of motion for the electronic dynamics. For this we turn to the 
Liouville-von Neumann equation for the evolution of the electronic density matrix 

ihp= [H e i,p] + ~1lp (14) 
n 

where the first term represents the unitary evolution of the uncoupled electronic states in the diabatic representation 
and the second term the non-unitary, dissipative dissipative dynamics due to the coupling to the phonon bath. Since 
the electronic state-space consists of roughly 100 energy levels, we restrict our attention to the population dynamics 
and decouple populations from coherences according to the Bloch model, 

Ti-aabb = ~ k a b + 5 a b ^ k ac (15) 

c 

Tlabab = r ^2( k ac ~ he) (16) 
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FIG. 1 Density and population of singlets (A) and triplets (B): density of states - black; initial population - blue; population 
at lOOps - red. 



where k a b are the rates of elementary interstate transitions. For internal conversions within the diabatic excited states, 
we assume that the phonons thermalize rapidly on the time-scale of the electronic dynamics such that the one-phonon 
transitions rates can be determined within the Markov approximation 



kgb = 7T 



(g n y) 2 

fru>e 



(71$ + l)(r(w£ - UJab) - T(tjj£ + UJab)) 



(17) 
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FIG. 2 Computed Condon spectral densities for (PPV) 32. Absorption (blue) and red-shifted emission (red) spectra between 
the adiabatic ground state and the excited states and the black curve is the oscillator strength for vertical excitation from 
the ground state. Spectral features correspond to excitation of C=C vibronic stretching modes in the excited state (upon 
absorption) and in the ground state (upon emission) as detailed in the inset figure. 
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FIG. 3 E-h capture population dynamics of essential states with and without e-h symmetry, for singlet (A) and triplet (B) 
cases. Solid curves are with e-h symmetry enforced and dashed curves are for broken e-h symmetry. In both figures, we show 
the population dynamics for the lowest two excitonic and charge-transfer states. 



Here nj is the Bose-Einstein distribution of phonons at T = 300K and T is the empirical broadening of O.OleV. Note 
that for an elementary a — > b population transfer to occur there must a phonon mode lu^ at that transition frequency 
Lu a b- Thus, phonon mediated dynamics is restricted to the excited state manifold where the energy level spacing is 
commensurate with the phonon energies. In contrast, single-photon processes occur across the S — > S± optical gap 
and the rate of spontaneous emissions is given by 



^ (1 + nK0) fe 



fc Q o = 7^(l + nKo))7Hf (18) 



where fi a o are the transition dipoles of the excited singlets. These we can comput directly from the Wannier functions 
or empirically from the photolumincscence decay rates for a given system. Photon mediated transitions between 
excited states are highly unlikely due to the density factor of the optical field. In essence, so long as the non- 
equilibrium vibrational dynamics is not a decisive factor, we can use these equations to trace the relaxation of an 
electronic photo- or charge-transfer excitation from its creation to its decay including photon outflow measured as 
luminescence. 



III. ELECTRON-HOLE CAPTURE KINETICS 

A. Capture cross-section and chain-length 

In modeling of c-h capture, propagation starts from the free e-h pair represented by the charge-transfer (CT) 
configuration with an electron and a hole on the far ends of the polymer chain, for example: 

|Vo) = (e - - ©) (19) 

We considered chains with n = 2,4, 8, 16, and 32 repeat units corresponding to longer and longer polymer chains. 
The salient data from each of these calculations is presented in Table |I] In Fig|T] we show the diabatic density of 
states (DOS), and the spectral distribution of the starting configuration in S and T state-spaces for a 32 repeat unit 
chain (PPV)s2- Singlet and triplet populations arc also given at an intermediate stage of relaxation at lOOps. In 
both cases, the system evolves to a metastable superposition state, where half of the density remains locked in CT 
states (Sct or Tct) above the lowest excitons (Sxt or Txt) - (2a) The branching into CT and XT channels is due to 
e-h symmetry, which separates excited states into even and odd representations under e-h transposition. Thus, XT 
states, which are even, are not vibronically coupled to CT states, which are odd. As shown in Fig|2K, Sxt and Txt 
are formed at different rates. Whereas Sxt, Sct and Tct show almost parallel population growths, formation of 
Txt is far slower. 

Assuming that after an initial approach time, the formation of the lowest energy excitons follows first-order kinetics 
corresponding to the decay of some precurser state, we can fit the population data to an exponential form and obtain 
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the formation rate constants. For the longest chain considered {{PPV)^), the formation half-times are t(Sxt) = 
17.80ps and t(T X t) = 174.77ps,i.e., r = 9.8 corresponding to a 77% EL/PL efficiency, which is systematically higher 
than observed in both LEDlJj; 0) and PADMRI0) data, but nevertheless consistent with the observed trend. For 
physical systems in which disorder limits the effective conjugation length, one might expect the effective conjugation 
length to be more on the order of 8 to f repeat units. For PPV this corresponds to r f=a 3 which is more in line with 
the experimental observations. 

The intrinsic distinction between S and T e-h captures is readily understood in terms of different exciton binding 
energies^ El EI 0) £b = £ct ~ £xt)- As seen in Fig^ the singlet excitonic band SjfTiiearly overlaps with the 
charge-transfer band Sct, whereas in the case of the triplet exciton, T\ lies at the bottom of a separate band of 
bound e-h pairs, about f .5eV below the CT continuum. Thus Txt formation requires on average a longer sequence 
of vibron-mediated transitions, each falling into the phonon frequency range. 



B. Breaking e/h symmetry 

The higher efficiency of S recombination becomes more apparent when e-h symmetry is lifted. We slightly broaden 
the conduction band and squeeze the valence band, so that /// = —1.1. The small e-h asymmetry changes negligibly 
the electronic spectrum, but opens a weak vibronic channel for CT— >XT internal conversions. The resulting population 
dynamics are shown in Fig[2j3- Here we see a fast build-up of low-lying XT and CT states in the first 100 to 200ps, 
followed by a slow conversion of CT states into excitons. The initial dynamics occurs approximately without e-h 
parity crossovers and the formation rates of low-lying XT and CT singlets and triplets are about the same as in the 
symmetric case. However, the triplet XT to CT branching ratio decreases drastically to about 1:2 and suppresses the 
formation of Txt- Moreover, subsequent Tct Txt relaxation is very slow due to the large binding energy of triplet 
excitons and the low density of states between Tct and Txt- Thus Txt population reaches only 40% after 800ps 
propagation. In contrast, Sxt is slightly favored over Sct in the initial capture, and further Sqt — * Sxt conversion 
occurs on a time scale of about 80Qps. l(20(l 

Note as well that Sct and Tct are very close in energy because of the lack of appreciable exchange between 
separated electrons and holes. Hence, intersystem S — > T crossing of long-lived CT states due to spin-orbit coupling 
(not included in the present model) is highly likely to occur prior to final e-h binding. 



C. Analytic Model-Universal Scaling Relation 

Let us consider the relaxation from the lowest unbound (CT) state to the lowest bound excitonic (XT) state due 
to coupling to the lattice vibrations. Since the exciton binding energy is greater than energy for a single phonon 
transition, we need to think about this a multi-phonon process. If we assume that this occurs as a set of discrete hops 
of energy SE « fuvphonon-, then the binding energy is eb — N5E. Furthermore, if we assume that the hopping rate 
between adjacent vibronic levels is the same for all steps, then the population of a given vibronic sublevel is 

hj = k(rii-i - rii) (20) 

where we assume that downward transitions are favored over upward transitions. For systems at low temperature, this 
will be the case so long as ksT < SE. Looking at the population of the state at the bottom, nxr, which corresponds 
to the final excitonic state, we can easily solve these coupled kinetic equations 

nx T (t) = l-e-^^. (21) 



TABLE I Collected data for various chain lengths (symmetric case) . Here, rs and tt are the formation half-times (r = ln(2) / k) , 
r = ks/kt, and <j> = r/(r + 3) is the electro-luminescent efficiency. 



n 


E s - Et (eV) 


TS(PS) 


r T (ps) 


r 




2 


1.552 


29.3541 


20.4905 


0.698 


0.19 


4 


1.309 


23.0852 


30.4756 


1.323 


0.31 


8 


1.199 


18.8897 


49.8341 


2.646 


0.47 


16 


1.174 


17.1152 


60.9348 


3.560 


0.54 


32 


1.146 


17.8005 


174.771 


9.818 


0.77 
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FIG. 4 A. Linear relation between polymer chain length, n and the rate-constant ratios for singlet and triplet cases. B. 
Singlet- Triplet gap vs. 1/n. Solid curve is fit to quadratic form. 



Now we ask, what is the time required for 1/2 of the initial population to appear in the lowest state starting from the 
highest? Not surprisingly, this scales with the number of intermediate states and consequently, the binding energy. 
Thus, 



which is consistent with our simulations discussed above. 

Next, taking the singlet-triplet energy gaps into consideration. To a lowest order approximation the singlet-triplet 
energy difference is approximately equal to twice the electron/hole correlation energy AE$t = IK. Using the data 
in Table 1, we can fit AE$t to a simple functional form 



where K x = 1.13eV and = 0.55eV and = Q.30eV, where n is the chain length. Thus, 



el 4 + AE 
r oc —2- = —2 — _ 



£ 



B 



Finally, using e§ oc 1/n, we obtain 



AE ST = Koo + K {1) jn + K {2) /2n 2 , 



el + K^+K^/n + K^/n 2 

' B 

r oc (1 + K {1) + nK^ + 0{K {2) /n). 



(23) 

(24) 
(25) 

(26) 



Since K^ is twice the correlation energy of an electron/hole pair in an infinitely long chain, this we can assume that 
for organic polymers, its value is more or less system independent. This universality consistent with the experimental 
observations reported in Ref.|@) Thus, we argue that the recombination process proceeds via a cascade through a 
series of intermediate states coupling the lowest charge-transfer states to the lowest excitonic state. 



IV. CONCLUSION 

In summary, we have simulated phonon-mediated intrachain e-h recombination in LCPs. The approach allows us 
to examine EL quantum yield in LCPs further case by case by combining first-principle description with available 
spectroscopic data. While we focus here on PPV and adjust el-ph coupling to its absorption spectrum, results can be 
generalized to most nondegenerate polymers. Apart from spin-degeneracy statistics, efficiency of singlet e-h capture 
outweighs the triplet one as a natural outcome of the higher binding energy of triplet excitons. For both S and T 
processes, the relaxation mechanism involves two steps as a result of approximate e-h symmetry. In the initial e-h 
capture, both low-lying excitons and same-spin CT intermediates are formed. The XT to CT branching ratio may 
vary, but is inevitably more favorable for singlet excitons than for triplet ones simply due to energetics. Relaxing 
to the higher lying singlet exciton requires fewer elementary relaxation steps than the lower lying triplet exciton. 
Subsequent conversion of CT states into excitons is forbidden by c-h symmetry and occurs on a much slower time 
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scale once this symmetry is broken. The weak e-h asymmetry of real conjugated systems favors again the singlet 
CT — ► XT relaxation by a factor, commensurate with the ratio of T to S exciton binding energies. 

What is desired is experimental detection of the transient kinetics we predict in out model. Ideally, the experiment 
would monitor the fluorescence signal following electron-hole injection onto a polymer chain. This is an extremely 
experiment since electro-injection typically entails multiple polymer chains with disordered morphology. This is a 
significant lacuna in our current model and we have begun work on incorporating multiple polymer chains into our 
theory to model inter-chain e/h transfer and recombination dynamics. 

What is also curious is that our predictions, a and the observation that the EL efficiency is dependent upon the rate 
of formation of emissive states, then the formation rate constant should be the limiting factor in the recombination 
process following electro-injection and not the electron/hole diffusion time. Surprisingly, the diffusion times for charge 
carriers is significantly longer than the recombination rates reported herein and for that matter predicted by other 
groups. Moreover, since carrier diffusion is more or less independent of spin-state (due to the fact the carriers are 
non-interacting), one expects that the EL efficiency to not depend upon the recombination process and should reflect 
the 1:3 ratio of singlets to triplets in the system. We hope that further experimental and theoretical work will shed 
some additional light onto the electroluminescence mechanism in organic polymers. 
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